Edge States at the Interface between Monolayer and Bilayer Graphene 
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The electronic property of monolayer-bilayer hybrid graphene with a zigzag interface is studied by 
both the Dirac equation and numerical calculation. There are two types of zigzag interface stacks. 
The dispersion and local density of states behave quit differently along the interface at the Fermi 
energy due to the different locations of the edge state. We hope our study can give some insights in 
the understanding of the transport and STM experiments. 
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I. INTRODUCTION 

In recent years, the experimental accessibility of the 
single and multi layered graphene sampleai"— has at- 
tracted considerable theoretical and experimental atten- 
tion due to its unusual electronic structure described 
by the Dirac equation, namely electrons in monolayer 
graphene have linear dispersion thus behave like mass- 
less Dirac fermions at the corners of the Brillouin zone 
(BZ)^. In the presence of magnetic field perpendicular to 
the graphene plane, the system shows anomalous integer 
quantum hall effector— (IQHE) which is different from 
that of the conventional two-dimensional electron system 
in semiconductor heterostructures. The Hall conductiv- 
ity in the IQHE of graphene shows half integer plateaus 
instead of integer ones in the IQHE of a normal semicon- 
ductor due to the different degeneracy at = Landau 
level. 

Edge states in graphene have been the focus of much 
theoretical study because of the important role they 
play in transport.— It is well known that there are 
two basic types of edges in graphene, namely, the arm- 
chair and zigzag edges. Some theoretical worki^r^ on 
the electronic structure of finite-sized systems, either as 
molecules or as one-dimensional systems, has shown that 
graphene with zigzag edge has localized edge states at 
the Fermi energy, but those with armchair edges do not 
have such state. Therefore, the zigzag edge states is es- 
sential to the transport properties in graphene since the 
localized edge state has contribution to the conductivity. 
However, most of the work were focused on the uniform 
monolayer and bilayer graphene. On the other hand the 
hybrid edge structure composed of partial monolayer and 
partial bilayer graphene , which is quit general in reality, 
has received not so much attention. Experimentally^^ 
the anomalous quantum oscillations in magnetoconduc- 
tance was observed due to the peculiar physics along the 
interface. Some of the transport properties and the pres- 
ence of interface Landau levels was explored within an 
effective-mass approximationi^2i2^ 

In this work2^, we study the electronic properties of 
the hybrid interface systematically via both tight bind- 



ing model and its effective theory in the continuum limit 
- the Dirac equation. The edge states in graphene can 
be studied experimentally by using a local probe such 
as scanning tunneling microscopy (STM). The STM ex- 
periments measure the differential conductance which is 
proportional to the density of states. Thus we study 
the local density of states (LDOS) in different hybridized 
zigzag edge graphene. It is defined as: 



N{r,eV) = \-^^{r)fS{eV-E^), 



(1) 



where ^'^(r) is the eigenfunction with energy The 
LDOS shows the strength of the local electronic density 
which is related to the strength of the signal in STM 
experimental data. Therefore, the LDOS can show us 
the signals of the different edge state as in our previ- 
ous worki^i In this paper, our study of the step edge 
graphene shows that there are two different zigzag step 
edge, and there always exists zero energy states localized 
near the zigzag step edge but the distribution of LDOS of 
these edge states strongly depends on the details of how 
the edge stacks together. We find the energy dispersion 
around the Dirac cones also presents different characters 
for different edge arrangements either within or without 
a magnetic field. 

This paper is arranged as follows: In section II, we set 
up the model Hamiltonian in different geometries. The 
zero energy solutions in zero field and in magnetic field 
are obtained by solving the Dirac equations and numer- 
ical diagonalization is implemented in a finite system in 
section HI and IV, and some discussions and conclusions 
are in the section V. 



II. MODEL AND DIFFERENT GEOMETRIES 



We consider the bilayer graphene with AB stacking, 
as shown in Fig. (TJ its tight-binding Hamiltonian can be 
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FIG. 1: (Color online)The schematic pictures for two kinds 
of monolayer-bilayer interface. The atoms of the extended 
bottom layer (layer 1) are indicated by black dots while the 
green dots represent the top layer which terminates at the 
interface. In the left plot, the top layer is ended with B2 sites 
il.e.b) while on the right it is ended with A2 sites (h.e.h). 



written as: 
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H t ^ ^ ^ ^ ^i;m,nibi-^m,n ~i~ ^z;m — l,n ~1~ ^2;m,n — l) 

i—1 m.n 

- t±J2"'2-m.,nKm,n+h.C., (2) 



the time reversal symmetry, the two Dirac points are not 
equal to each other, and both must be studied. 

According the lattice orientation we adopt which is 
shown in Fig. [TJ the zigzag interface is along the x direc- 
tion. For simplicity, we consider an infinite stripe along 
the X direction, therefore the system has translational 
symmetry along the x direction thus remains a good 
quantum number. We then do Fourier transformation in 
the X direction and reduces the 2D problem to ID. There 
are actually TWO distinct geometries which are physi- 
cally different, i) The outmost sites of the upper layer 
are the B2 sites which do not stack directly on the lower 
layer atoms as showed in Fig. [Ijleft). The B2 sites are 
the low energy degrees of freedom (along with Ai ) which 
are kept if one further considers an 2 x 2 effective theory 
on energy scale e ^t±. We label it the low energy sites 
boundary {l.e.b.); ii) the A2 sites are the outmost sites 
on the upper layer as showed in Fig. [IJright). The A2 
sites, together with the Bi sites which they stack right 
on top of, form the dimer sites. In the 2x2 low energy ef- 
fective theory the wavefunction have almost zero weight 
on those dimes sites when e <c As a result these 
dimers are ignored in this limit. In other words, they are 
occupied considerably only at high energy (comparing to 
t±). Therefore, we shall address such interface as the 
high energy sites boundary [h.e.h.). 



where ai-m,n {bi-m,n) is the annihilation operator at posi- 
tion (to, n) in sublattice Ai (Bi), and i = 1, 2, indicating 
the two layers. The first term is the Hamiltonian within 
each layer, and the second term describes the interlayer 
coupling in which we only consider the hopping between 
the two atoms stacked right on top of each other. Let 
us label the bottom (extended) layer as layer 1, the half 
plane upper layer as layer 2. In this work we will only 
consider the A — B (Bernal) stacking. One can expand 
the effective Hamiltonian near the two Dirac points K 
and K' which are time reversal symmetric partners. In 
momentum space, the hamiltonian near K can be written 
as: 



(3) 



where 



Hk = 



in which k 

{ai;k,bl-k,0'2;k 

before, Hk' - 




b2-/)- 
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t_i_/vF, 

For the other Dirac point, as stated 
Here we only consider the zigzag 
type interface(or edge) to explore the localized edge state. 
Without magnetic field and the interface, it is sufficient 
to discuss just one Dirac cone in the continuum model 
due to the symmetry. But with the interface breaking 
the inversion symmetry and the magnetic field breaking 



III. INTERFACE PROPERTIES IN ZERO FIELD 

For a semi- infinite sheet, it is well known that the ex- 
istence of the zero energy edge modes in both mono- 
layer and bilayer graphene. Presumably, such modes 
are also expected at the interface between them. Let us 
consider the following geometry: a half plane of mono- 
layer graphene and a half plane of bilayer graphene joined 
along the zigzag edge; and look for solution(s) with zero 
eigenenergy by using the Dirac equation. 

Before we present the results, let us discuss the bound- 
ary conditions at the interface firstly as in Ref. [sil Let 



and 



*bi(a;,2/) 



tpB{x,y) 



ipA2ix,y) 
\i^B2{x,y)j 



(5) 



(6) 



be the wavefunctions at Dirac point in the monolayer 
and bilayer respectively. Suppose the interface locate at 
?/ = 0, the boundary condition for monolayer is then 
straightforward: both components of the wavefunction 
must be continuous. For the upper layer, it terminates 
at y = and therefore satisfies the open boundary con- 
dition. Note that the term 'terminates' indicates the last 
row of lattice sites are the high/low energy sites, how- 
ever, the boundary condition is not the wavefunctions 
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being zero on these sites. They should be the wavefunc- 
tions on sites one unit ceU 'outside' the boundary being 
zero. To sum up, we have: 



Ua2{x,0)^0 l.e.b. 
\i'B2{x,0)^0 h.e.b. 



(7) 



In the monolayer region with zigzag interface, we do 
the substitution k - 
k 
k* 



kx + dy hi the Dirac hamiltonian 
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The zero energy solution is 



(t>B{y) 



(8) 



In analogy, we do the same substitution in the bilayer 
Hamiltonian (Eq. S]) and obtain its zero energy solution: 



t>Biiy) 

i>Ai{y) 
i>A2{y) 



(9) 



Applying the requirement that the wavefunctions re- 
main finite at infinity, for the l.e.h. interface (so 4'A2{y = 
0) =0), one easily find that nonzero solutions only exist 
for kx > Q- The solution is 

Ci = C2 = Ai = A:i = Ai = 0, A2 = const. (10) 

For the h.e.h. interface, one finds that for > Q 

Ci = C2 = Ai ^ A2 = A3 ^ 0, A4 = const. (11) 

For the other Dirac cone, the solutions remain the same 
but only exist for kx < 0- The only nonzero constant is 
to be determined by normalization. Both results are in 
agreement with tight-binding analysis^. 

Even though the zero energy states exist for both types 
of interface, there is an important difference between 
them. In the l.e.b. case, the wavcfunction is only none- 
zero in the upper layer, which is trivial as such mode is 
expected when a graphene sheet is terminated at a zigzag 
edge. However, the h.e.b. is less trivial. The wavefunction 
also lives on the extended layer at the interface where no 
cut is present. We interpret this as, when the dimer sites 
arc the boundary, the interlayer coupling t± imposes an 
energy cost for electrons to go through the interface in 
the extended layer which can be considered as an effec- 
tive potential barrier. The potential barrier can localize 
the electron states along the interface. 

We numerically diagonalize a system with a finite 
width up to 600 unit cells in the y direction. The in- 
tralayer hopping strength t is set to identity and the 
interlayer hopping strength t±_ = 0.2t. The dispersion 
relation is shown in Fig. [21 Compare the dispersions in 
different geometries, the h.e.b. edge has an obviously 
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FIG. 2: The dispersion around Dirac cone for the step 
graphene with (a)l.e.b and (b)h.e.b respectively. 



stronger level anti-crossing feature than the other. The 
reason is that for l.e.b. edge, the zero energy edge state 
just locates on the bilayer graphene which has quadratic 
dispersion and has nothing to do with the monolayer part 
with a linear dispersion. However, in the case of h.e.b., 
the zero energy edge state also has component on the 
monolayer graphene at the interface, the linear dispersion 
part for monolayer should also connect to the edge state 
which induces more energy level anticrossing around the 
Dirac points. Figl3] shows the LDOS for the two kinds 
of step graphene. We label the coorindinate of the ex- 
tended layer by d € [0, 600) and the d e [900, 1200) for 
upper layer. Therefore, the interface locates at d = 300 
and d = 900 for layer 1 and 2 respectively. One can no- 
tice that the localized edge state shows as a peak at zero 
energy along the interface. The prominent difference is 
that the peak just appears on the top layer of bilayer part 
in the l.e.b. case; but for the h.e.b. case, the zero energy 
peak appears on both two layers although still locates 
at the bilayer side. The numerical results are in agree- 
ment with the analysis by the Dirac equations and we 
conclude that the two kinds of edge arrangements should 
have different consequences in experiments since the dif- 
ferent distributions of the zero mode. Here we notice that 
the other peaks in LDOS at the end of the finite system 
is the signal of the general monolayer or bilayer zigzag 
edge graphene as discussed in many others woihr^''— 



IV. INTERFACE PROPERTIES IN MAGNETIC 
FIELDS 

In the presence of magnetic field, by minimal cou- 
pling the Dirac equation should be modified by doing 
the substitution k k + Assume the magnetic 

field is along the direction of perpendicular to the plane 
B = Bz, B > 0. We adopt the Landau gauge here 
A = (Ax, Ay) = {—yB, 0) due to the translational invari- 
ant along X direction. Thus the Hamiltonian of mono- 
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FIG. 3: (color online) The LDOS of the step bilayer graphene 
with l.e.b.(left) and h.e.b. (right) respectively. The width of 
of monolayer L = 600. In plot, the region d £ [0 : 600) stands 
for the lower extended layer and d £ [600, 1200) is the upper 
layer. 



layer becomes: 



c 



eBy 



and its zero energy solution is 



My) = C2e-''-y+'^y" 



(12) 



(13) 



Similarly, the zero energy solution for bilayer graphene 
becomes: 



My) ^ A.e-'^^y+j^f 
(l)A2{y) = A-ie^-y-'&y' 

AB2{y) = {-Ai-iy + Ai)e-''-y+^y" 



(14) 



Apply the same boundary conditions as in the zero field 
case, we will have the solutions. For Le.6., one gets 



C2 = ^1 = ^3 = ^4 = 0, 

For h.e.h. the solution is 



Ci = = const. (15) 



C2 = Ai = A4 = 0, Ci = A2 = consti, A3 = const2. 

(16) 

The constants are to be determined by normalization 
conditions. One immediately notices that the wavefunc- 
tion lives only on the A sublatticc. We should note that 
on the other Dirac point, the solutions remain the same 
form, but resides on the B sublatticc. Another impor- 
tant feature of the solution is that for h.e.b. we actually 
have TWO independent solutions here. 



A. Dispersion Relation 

The general problem in magnetic field can be written 



as 



(17) 



for the monolayer while that for the bilayer is 





d^+^ 

J_ I -9^ + ^ 7 

V2 7 d^+C 

+ C 



hi 



where ^ 

7 



VF^c 

the monolayer is 



(18) 

It is known that in the bulk, the solution to 



eB 



4-. 



(19) 



where the e = ±VN, TV = 0, 1, 2 . . . , and £> ^'s are the 
parabolic cylinder functions, which combined with the 
factor l/r(u + 1) give us the eigen wavefunctions of a 
harmonic oscillator ipi^{£,). The bulk solution to the bi- 
layer Hamiltonian can be written in a similar way: 



£M2±l)w,. 



(20) 



where e = ± V i+^-+^^-±V(i+y)-+! l!2 (j- ^ 0, 1, 2, . . . ) is 
the eigen energy. 

However, for the interface problem, we need solutions 
on the half-plane. In this case, we let the bilayer live on 
the y > side, so the monolayer is on the y < part. 
Therefore, for the bilayers, the solutions on (0, 00) take 
on the same form as in the bulk, but fs are no longer 
required to be integers. Instead, we now should replace 
j's by 



1 



1 



(21) 



and e now varies continuously. For the monolayer, the 
solution on (— oo,0) can be chosen as 



(22) 



With the above solutions, and combined with the bound- 
ary conditions that are already discussed in the zero field 
cases, we obtain the transcendental equations that dic- 
tates the dispersion relations. For the l.e.b. interface, 

(e^ - (ji + l))Dj,^i{~V2k^){D,2{V2k.,)D,,{-V2k^) 

+D,2_i{V2k,)Dj,+i{-V2k,)) 
= (e' - (.72 + l))Dj,^i{-V2k,){D,2iV2k,)D,,{~V2k,) 

+D,2_,{V2k,)D,,+,{-V2k,)), (23) 
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where we have set = 1. For the h.e.b., the dispersion 
equation is 

(e' - (.71 + 1) - j^)D^,{~V2k,){D,2{V2k,)Dj,{~V2k,) 

+D,2_i{V2k^)D^,+i{-V2k^)) 
= (e' - in + 1) - -/^)D,,i-V2k,){D,2{V2k,)D,,{~V2k,) 

+D,2_i{V2k,)Dj,+i{~V2k^)) (24) 

However, one must note that in the presence of mag- 
netic field, the time reversal symmetry is broken, there- 
fore, the other Dirac cone, the time reversal partner, no 
long behaves the same way. So the dispersion relation 
must be calculated separately. By the same approach, 
one can get, for the l.e.b., 

- n)Dj,{^V2k^){j2D,2{V2k.,)Dj,_,{~V2k^) 
+e^D,2_^{V2k,)D,,{~V2k,)) 
= {e^- - j2)Dj,{-V2k^){jiD,2{V2k^)Dj,-ii-V2k^) 
+e^D,2_,{V2k.,)D^,{~V2k,)); (25) 

for the h.e.b., 

(J2 + l)(e' - Ji - i'')D,,{-V2k^){32D,2{^k^) 

D,,_,{-V2k,) + e^D,2_^{V2k,)D,,{^V2k,)) 
= in + l)(e2 - j2 - ^^)Dj,i~V2k^){hD,2{V2k^) 
D,,_ii-V2k,) + e^D,2_,{V2k,)Dj,i-V2k,)). (26) 

It is easy to see from the dispersion equation that when 
\kx\is very large, the eigen-energy should restore to either 
the monolayer value or the bilayer value as the parabolic 
cylinder functions only converge to zero for integer 
v at infinity. But what interests us is how these bulk 
Landau levels connect with each other when crossing the 
interface. To study that, we solve the above equations for 
e and k^ near the interface fc^, = for a finite size sys- 
tem as shown in Fig. S) We see that at the two Dirac 
cones, the Landau levels match in different ways. Fur- 
thermore, we present the numerical result for directly 
diagonalizing a ribbon with width L = 600. Our diag- 
onalization results consistent with the analytic analysis 
also. Fig. [5] and Fig. IH] shows the energy spectrum and 
LDOS for step bilayer graphene in a magnetic field with 
l.e.b. and h.e.b. respectively. In magnetic field the en- 
ergy spectrum is split into a set of Landau levcls(LLs) 
and the Landau level splitting in bilayer is smaller than 
that of monolayer graphene since the energy of LLs in 
bilayer system satisfies E oc n{n — 1) and E oc ^/n in 
monolayer. Therefore, LLs are not matched at the inter- 
face. 

The major difference between the two sets of disper- 
sion relation is that the Landau levels in the two distinct 
regions connect in two different ways. The dispersion 
at the interface with l.e.b. is more flat than that with 
h.e.b.. These interface Landau levels was discussed by 
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FIG. 4: (color online) The dispersion around Dirac cones in 
the presence of magnetic field. The upper two figures are 
for the l.e.b. and h.e.b. interfaces around one Dirac cone 
respectively, and the lower two figures are that for the other 
Dirac cone. 



M. Koshino et.al. within an effective-mass approxima- 
tion2^. From the plot of LDOS, the same as in the case 
of zero field, according to the different distribution of the 
edge state, the presence of the peak a.t E = along the 
interface only on the upper layer for l.e.b. and on both 
two layers for h.e.b.. 

We also use the Dirac equation result to study how 
the dispersive Landau levels near the interface changes 
when the magnetic field varies. If we measure the energy 
in unit of hiOc, the monolayer bulk Landau levels remain 
unchanged as the magnetic field is being tuned; but for 
the bilayer, the effective interlayer coupling 7 ^ l/wc so 
the bulk Landau levels vary with the magnetic field. In 
strong field 7 ~ 0, so the two layers behave as if they 
were decoupled, and the Landau levels' energy becomes 
just like the monolayer's but with a two-fold degener- 
acy. I n weak fi eld limit, 7 ~ 00, in this case one gets 
e ~ \/n{n — l)/7. The dispersions for 7 G 0.2,1,2,5 of 
both l.e.b. and h.e.b. at one Dirac point are shown in 
Fig. [7] and Fig. [5] respectively. Here we only plot the 
Landau levels up to n = 6. The evolution of dispersion 
on the other Dirac point have the similar behavior ex- 
cept different energy level connection as shown in Fig. |4l 
As the strength of the magnetic field is increased, the 
bilayer LLs become more denser which induces more me- 
diate states around the interface. The dispersion of the 
l.e.b. interface is always more flat than that of h.e.b. in- 
terface while tuning the magnetic field. 
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FIG. 5: The energy spectrum and LDOS of the step bilayer 
graphene with l.e.b. in a magnetic field. The strength of the 
magnetic field is expressed as magnetic flux (j> in each unit 
cell. Here we set the magnetic flux per each unit cell (f) — 
(?!>o/1315 ((^0 is the magnetic fiux quanta) which corresponds 
to B ~ GOT. 



V. DISCUSSION AND CONCLUSION 

In this paper, we discuss the edge state at the inter- 
face between monolayer and bilayer graphene. Physi- 
cally, there are two types of interface structures named 
by the l.e.b. and h.e.b. due to the different terminations 
at the interface. By studying the system with both effec- 
tive theory and numerical calculation, we find that with 
or without magnetic field, zero energy edge states exist 
for both types of interface, however, the LDOS shows 
different features. For l.e.b. the zero energy edge states 
only live in the upper layer which terminates at the in- 
terface while the h.e.b. induces an enhanced LDOS in 
the bottom layer. Another major difference between the 
two different geometries is that the dispersion of h.e.b. 
system shows a stronger anticrossing feature in the zero 
field. Both differences can be interpreted as the follow- 
ing. When an electron goes through the h.e.b. interface, 
it can choose to lower its energy through the interlayer 
coupling. This can be viewed as the h.e.b. imposes an 
effective potential barrier at interface in the extended 
layer. Since the localized zero energy edge states locates 
on both two layers, the monolayer dispersion which is 
linear should smoothly connect to the edge state in the 
h.e.b. case. This can explain the strong energy level an- 
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FIG. 6: The energy spectrum and LDOS of the step bilayer 
graphene with h.e.b. within magnetic field. The parameters 
are the same as in Fig. \5\ 
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FIG. 7: The dispersion around Dirac cones K of l.e.b. type 
interface in different magnetic field with 7 = 0.2, 1,2,5 re- 
spectively. 
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7=0.2 7=1 
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FIG. 8: The dispersion around Dirac cones K of h.e.b. type 
interface in different magnetic field with 7 = 0.2,1,2,5 re- 
spectively. 



pcrsion of Landau levels continucously goes through the 
interface in different manners at the two Dirac cones. 
At Dirac cone K{K'), the zero energy Landau levels in 
the bilayer region split, one branch rises up and becomes 
the n = 1 Landau levels in the monolayer while the 
other branch remains as zero energy states, and the other 
higher Landau levels continues from bilayer to monolayer 
accordingly, i.e. n = 1 , . — > n = 2 .... However, 

^ I bi I mono ' 

near the other Dirac point K'{K), the zero energy state 
remains intact, so the other higher Landau levels connect 
as n = l|, . — s- ri = l| .... We also sec an effect of 

I bi I mono 

the effective potential barrier imposed by the h.e.b. in 
the dispersive Landau levels near the interface. For a 
certain range of field strength the dispersion develops a 
local maximum for h.e.b. while the l.e.b. only develops a 
plateau feature. 

In summary, the physical properties of hybrid graphene 
systems are mostly dominated by both that of mono- 
layer and bilayer. The Landau levels remain the same 
away from the interface. But the existence of dispersive 
Landau levels near the interface could be related to the 
unexpected feature other than that of the monolayer and 
bilayer graphene in the magneto transport experimenli^^. 
Further study of the hybrid structures for a more realis- 
tic setup, like including edge disorder, gate voltage, etc., 
are needed to undertand the experimental data. 



ticrossing in the h.e.b. interface. Similar differences be- 
tween l.e.b. and h.e.b. interfaces are also discussed for 
transimission coefficients across the interface^'^'^. For 
the h.e.b. the transmission probability is reduced signifi- 
cantly for incoming electrons with energy E ^ t± which 
can be interpreted in a similar way. 

We also see that in presence of magnetic field the dis- 
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